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ABSTRACT 


Air  data  out;  Us  obtainable  from  pressure  measurements  on  a  h  -mispherical 
probe  have  been  im  ssligaled  analyti  ally  through  tin  Hypersonic  Mach  number 
range.  Angles  >f  at’  ck  from  +50  *  i.>  - 2.0 *  and  angles  of  sideslip  to  ±15°  arc 
considered  using  a  fiv.  orifice  probe  (one  centerline,  two  each  in  angle  of 
.  ttack  and  angle  of  sidcsi'p  plane-  Emphasis  is  put  on  the  hypersonic  regime, 
wherein  air  data  outputs  are  shown  to  be  obtainable  by  using  a  simplified  set 
of  equations.  Specifically,  vehicle  attitude  can  be  obtained  from  pressure 
inputs  alone;  true  air  speed  ca  •  be  ubtained  with  the  additional  input  of  free 
stream  density.  Pressure  expressions  used  to  obtain  attitude  angle  for  the 
hypersonic  regime  arc  found  to  be  acceptable  for  Mach  numbers  as  low  as  .  5. 
Expressions  for  determining  the  uncertainties  in  the  air  data  outputs  resulting 
from  both  pressure  measurement  erroi  and  simplifying  assumptions  used  in 
deriving  the  air  data  equations  are  presented.  Discussions  on  the  capability 
of  determining  vehicle  altitude  and  Mach  number  with  inputs  of  hemisphere 
pressure  and  free  stream  density  arc  given.  Finally,  air  data  output  errors 
resulting  from  changes  in  nose  shape  due  to  material  ablation  are  also  considered. 
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SECTION  I 


AIN  i.  AW iJUIj  11  UUN 


Air  data,  as  used  in  this  report,  means  information  describing  the 
altitude,  air  speed,  Mach  number  and  altitude  of  a  flight  vehicle.  In  this 
stodv,  assumed  air  data  inputs  are  pressure  mas suremenf 3  on  a  hemispherical 
probe,  with  and  without  the  additional  input  ol  free  stream  density.  The  probe 
in  r.  typical  fiignt  usage  could  well  be  the  nose  (on  the  order  of  a  half-foot 
diameter)  of  the  flight  vehicle,  and  it  will  be  assumed  to  be  free  of  pressure 
perturbation  due  to  other  vehicle  components. 

Instrumentation  and  techniques  used  to  obtain  air  data  at  subsonic, 
transonic,  and  low  supersonic  Mach  numbers  are  highly  developed,  a 
bibliography  cf  this  work  is  presented  in  Reference  (1).  At  these  speeds,  the 
conventional  appioach  has  been  to  equip  the  flight  vehicle  with  external  probes 
provided  with  sensors  which  measure  stagnation  pressure  and  indicated  static 
pressures,  local  flow  angles  and  total  temperature.  That  information  can  be 
related  to  free  stream  values  of  static  pressure,  angle  of  attack  and  sideslip 
and  temperature  to  provide  data  on  the  relationship  of  the  aircraft  to  the 
atmosphere  through  which  it  is  flying  as  well  as  data  on  the  atmosphere  itself. 

It  is  normal  practice  to  design  probes  and  locate  static  sources  on  the  vehicle 
which  by  themselves  introduce  minimum  errors  in  the  measurement  of 
stagnation  and  static  pressures  over  the  flight  range  of  interest.  Once 
installed  on  a  vehicle,  the  air  data  computer  is  designed  to  correct  for  the  so- 
called  position  errors  of  the  flight  vehicle. 

To  speeds  of  about  Mach  2  or  Mach  3,  these  errors  may  be  regarded  as 
perturbations  to  the  basic  measurement  and  corrections  as  determined  from 
full-  scale  flight  tests  and  wind  tunnel  tests  are  applied.  At  high  supersonic - 
hypersonic  Mach  numbers,  however,  the  conditions  arc  vastly  different. 
Hypersonic  flows  are  characterized  by  strong  shock  wave  systems  lying  close 
to  tin  body  surface;  hence,  air  data,  static  pressure,  for  exanq  "e,  cannot  be 
obtained  through  small  corrections  to  some  measurement.  To.al  or  stagnation 
temperature  measurement  is  subject  to  real  gas  effects  and  proper  inter¬ 
pretation  of  stagnation  enthalpy  from  this  measurement  is  difficult  (see 
Reference  2), 

In  the  past  half  dozen  years  or  so,  much  valuable  experimental  information 
has  been  obtained  concerning  air  data  at  supersonic  Mach  numbers  --  particularly 
with  regard  to  vehicle  attitude  on  X-15  flights,  see,  for  example.  References 
(3-6).  Attitude  angle  for  these  flights  was  measured  using  rotating  hemispherical 
null  seeking  pressure  probe.  This  investigation  was  undertaken  in  order  to 
describe  fully  the  maximum  utility  of  a  simple  fixed  position  pressure- 
instrumented  hemisphere  probe  at  high  Mach  numbers,  and  also  to  consider 
the  maximum  obtainable  air  data  outputs  obtainable  from  this  sensor  at  the 
lower  Mach  numbers. 

Although  the  study  considered  air  data  investigation  over  the  range  oi  Mach 
numbers  up  to  20  for  altitudes  to  30G,  000  feet,  because  of  present  air  data  needs 
for  hypersonic  cruise  as  well  as  lifting  re-entry  vehicles,  emphasis  was  put 
on  (he  high  Mach  number  range.  Furthermore,  of  the  air  data  outputs  defined 
above,  primary  emphasis,  again  because  of  present  needs,  was  given  to  vehicle 


attitude  (angles  of  attack  and  sideslip)  and  velocity  (true  air  speed).  Specifically, 
ttie  study  consisted  of  defining  the  e  uations  needed  to  describe  the  pressure 
distribution  over  a  hemisphere,  pa rticuj.;  i  !y  in  the  hypersonic  range,  and  then 
using  these  equations  developing  expressions  for  the  air  data  outputs,  (angle 
of  attack,  true  air  speed,  etc.).  Throughout  the  study,  frequent  use  of  error 
analysis  was  employed  in  order  to  assess  no!  only  the  usefulness  of  the  various 
expressions,  but  as  a  useful  tool  in  determining  orifice  location,  and  the 
validity  of  simplifying  assumptions. 

The  air  data  investigation  presented  herein  was  written  such  that  similar 
air  data  outputs  (c.  g.  oL  and  )  are  grouped  in  individual  sections.  Angle 
of  attack  and  angle  of  sideslip  are  discussed  in  Section  II,  True  air  speed 
and  indicated  air  speed  investigations  are  given  in  Section  III.  Mach  number 
and  altitude  because  of  their  similar  dependence  on  the  input  of  free  stream 
density  are  both  considered  in  Section  IV.  A  cursory  examination  of  the 
degree  of  nose  ablation  for  a  typical  flight  case  and  the  effect  of  this  ablution 
on  the  air  data  attitude  equations  is  giver,  in  Section  V.  Each  section  was 
written  with  the  intent  of  making  it  a  fairly  complete  discussion  in  itself  of 
the  particular  air  data  output  being  considered. 

Results  arc  given  in  forms  which  can  be  easily  assessed  for  accuracy 
and/or  modifications  for  a  particular  application.  Suggestions  for  additional 
work,  primarily  experimental,  to  further  develop  the  usefulness  of  this  air 
data  probe  are  given. 
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SECTION  II 


VEHICLE  ATTITUDE  ANGLE 


1. 


DEFINITION  OF  VEHICLE  ATTITUDE  ANGLE  AND  AXIS  SYSTEM 


Consider  n  body  at  some  arbitrary  attitude  with  respect  to  the  free  stream 
velocity  vector  Uw  ,  as  shown  in  Figure  1.  To  describe  the  attitude  of  the 
body  with  respect  to  this  velocity  vector  an  axis  reference  system  is  defined. 
This  system,  employs  the  body  axes  (see  c.  g.  Reference  3)  and  is  a  right- 
handed  orthogonal  act  with  the  origin  at  the  body  center  of  gravity  (r.  g. ).  A 
longitudinal  plane  of  symmetry  for  the  body  is  assu.ncd.  Then  the  axes  are: 

x  -  in  the  plane  of  symmetry,  directed  longitudinally'  forward. 

y  -  normal  to  the  plane  of  symmetry,  directed  along  the  right  wing. 

z  -  in  the  plane  of  symmetry  directed  "down.  " 

Further,  let  the  projections  of  the  velocity  vector  U  on  the  body  axes  be 
called  TT,  V,  w  for  components  along  the  x,  y,  and  z  axes,  respectively. 

The  angles  are  now  defined  in  terms  of  the  velocity  vector  and  its 
projections  on  the  orthogonal  body  axes.  /  lglo  of  attack,  OC  ,  anc  .ngle  of 
ideslip,  ^3  ,  are  defined  in  the  x-z  and  x-y  planes,  respectively,  as: 


.  -I  TT  -l 

angle  of  attack,  OC  *  cos  'fj  *  sm 


(ccSuc*) 


angle  of  sideslip, 
A  third  angle,  S 


Z3 


*  cos 


-1  TI 


(a'TuZpi 


-l  Tr 


(u 


Sift 


(w*  +  v*)v'*- 


(herein  unnamed)  that  is  useful  to  define  is: 

It 


s  = 


-l  (u.W») 
cos  — — - 


=  si  n 


-1  m 


u. 


(1) 

(2) 

(3) 


Th 

example. 


angle  S  has  sometimes  bear,  called  the  angle  of.  sideslip  (see,  for 
References  9-11).  However,  in  this  study,  the  sideslip  angle  is  the 
angle  defined  by  Equation  (2).  Note  that  while  Ot-  and  fb  are  defined  in  the 
planes  of  the  body  axes,  3  is  not.  Further  note  that  3  is  independent  of 
angle  of  attack,  i.  e.  the  vector  (u*  +  w z ) remains  constant  as  Of.  varies; 
whereas  is  net  independent  of  tX  ,  since  the  vector,  u,  varies  as  Of. 

varies.  For  this  reason,  <5  will  be  used  in  tire  derivations  of  an  expression 
for  CL  in  terms  of  the  pressure  distribution  in  the  x-z  plane,  wherein,  the 
pressure  distribution  in  the  x-y  plane  will  subsequently  be  shown  to  be 
independent  of  <6 

2.  ATTITUDE  ANGLE  DETERMINATION  AT  HYPERSONIC  MACH  NUMBERS 

The  capability  (including  accuracy,  range  of  use,  etc.  )  of  the  hemispherical 
pressure  probe  as  used  to  measure  angle  of  attack  and  angle  of  sideslip  depends 
on  the  following  factors: 

1)  The  accuracy  to  which  the  pressures  arc  known. 

2)  The  correctness  of  the  assumed  pressure  distributions. 


3)  The  choice  and  degree  of  exactness  of  the  derived  attitude  equations 
- -  tv  tn<j  » election  oj  oniicc  locations. 

Pressure  accuracvf,  of  course,  is  applicable  to  the  entire  Mach  number  range. 

nmnbc,  ro^°  faclor**  hnwevBr*  require  different  analysis  for  different  Mach 

on  the  huu  gr«\GS‘  <IC’  a3.statcd  ln  the  introduction,  emphasis  wiU  be  placed 
n  the  hypersomc  Mach  number  regime,  consideration  is  first  given  to  this 
range  and  then  to  the  lower  Mach  numbers.  K 

a.  Pressure  Accuracy 

i  .  ,  ,Thc  air  error  equations  developed  are  general  enough  to  allow' 

introduction  of  pressure  error,  of  various  magnitude,  however,  plus  or  minus 
ore  percent  pressure  errors  are  used  in  this  study.  One  percent  errors  arc 
bu  I?  represent  present  state  of  the  art  pressure  transducer  capability 
but  rather  are  used  solely  for  purposes  of  illustration  and  comparison.  Present 
systems  are  better  than  one  percent,  perhaps  one-fourth  of  one  percent,  thus 
conclusions  as  to  the  actual  usefulness  of  a  result  are  conservative. 

b.  Pressure  Distribution  at  Hypersonic  Mach  Numbers 

Pressure  data  obtained  on  a  hemispherical  nose  can  be  thought  of  as 
two  separate  data  inputs:  pressure  distribution  and  absolute  pressure  magnitude 
Drtornn-at.o,,  of  atfuce  angle,  «.  and  A  will  ho  to  be  independent  of 

absolute  magnitudes,  dependent  on  distribution  only. 

The  most  widely  accepted  approximate  method  used  to  predict  the 
pressure  distribution  or.  a  hemisphere  in  hypersonic  flow  js  the  simple  Ncw- 
tonian  impact  flow  theory.  It  is  useful  to  consider  the  assumptions  upon  which 
lh  .  theoiy  is  based  and  to  include  herein  its  derivation. 

Consider  the  streamline  prirni'ntprimi  tkr,  i  . 

_..-T.  /->  /  .  ,  .  - ***e>  - ..a  duuttp  at 

t  bei0Wi*  No  mcntion  is  made  of  any  shock  process 

associated  with  die  flow-;  it  is  assumed  that  none  exists. 


( 


Thl°r^  as  s umption  that  is  made,  and  this  constitutes  the  basis  of  the  theory 
is  that  tne  streamline  must  turn  an  angle  (90-0)  up-n  impact  with  the  body 

f?fnink8°ad0inSe'  al1  ^  momentum  normal  to  the  surface,  m  \J-,n  ,  is  lost 
to  the  body,  whereas  the  tangential  component  of  momentum  is  entirely 


conserved.  Therefore,  the  normal  force  on  the  surface  is 

^  =  dt  i'™'  5T 

whore  the  mass  fJow  rate  normal  to  the  surface  is 


dLt 


nt 


=/>„ 


U„,„  5 


for  area  S. 


(4) 


Notice  the  assumption  of  an  incompressible  fluid,  jO <m  =  Constant,  is  implied 
here. 


Thus, 

but 

is 


v.  acre 


vL  s 


n  *  p* 

c  U*  cos  0 

P  -  m/s 


so  the  pressure  at  0  , 


(5) 

(6) 


P  */3LUj  cos  6 


%  m~z  /°' 


t 


cos 


6 


(7) 


U 


to  CO 


*  is  the  pressure  acting  on  the  surface  at  0  .  This  is  the  Newtonian  flow 

prediction.  As  the  impact  angle  goes  to  zero, cos  6 - *1,  the  fluid  at  this 

location  stagnates  and  the  Newtonian  theory  simply  gives: 


K/<x  -  2  <s> 

'  0 

or  P  ,  the  stagnation  pressure,  is  equal  to  two  times  the  dynamic  pressure,  q. 
T1  alue  of  P  is  known  experimentally  to  bo  incorrect.  Quite  often  a 
moe  ed  Newtonian  Theory,  whereby  the  value  of  P/cos*  0  is  set  equal  to  a 
knt!  stagnation  pressure  P  at  &  3  0  ,  is  used  in  the  pressure  distribution 
prediction;  i.  e. ,  ° 


P  -  Ps  cos1  6  (9) 

The  vc  lity  of  this  pressure  distribution  equation  at  high  Mach  numbers  is 
next  considered  in  view  of  analytical  and  experimental  results. 

Experimental  pressure  data  taken  in  the  Cornell  Aeronautical  Tab- 
oratory  hypersonic  shock  tunnel  (Reference  12)  on  a  12-inch  diameter  hemis¬ 
phere  cylinder  at  =  14,000  fps  is  shown  in  Figure  2.  These  results  are 
averaged  data  for  approximately  14  tunnel  runs,  and  are  also  shown  corrected 
for  the  effects  of  the  conical  flow  field  in  which  the  model  was  tested  using  the 


*  Sut  fe sated  by  L.  Lees,  IAS  Preprint  JJ4,  1955. 
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characteristic  conical  flow  solutions  of  Reference  (13),  In  conjunction  with  the 
experimental  tests,  a  computer  pro^iam  for  a  real  -  reacting  gas  was  used  to 
predict  the  hemisphere  pressure  distribution  for  the  test  conditions  of  the 
program.  Additionally,  exact  numerical  solutions  for  an  id  :al  gas  (Reference 
14)  arc  shown  in  Figure  2. 

Several  important  points  are  evident  from  Figure  2: 

1)  The  corrected  experimental  data  and  the  real  gas  solution 
at  Ujo  =  14,  000  fps  arc  both,  in  excellent  agreement  with  the 
ideal  gas  solutions. 

2)  The  ideal  gas  solutions  are  practically  invariant  with  Mach 
number  for  6  <  M  <  30. 

•y 

3)  These  curves  lie  below  the  cos'"  Q  curve,  showing  an 
inadequacy  in  using  the  modified  Newtonian  distribution. 

In  view  of  these  results,  the  numerical  real  and  ideal  solutions  and  the 
experimental  data  are  accepted  as  presenting  the  correct  hvpersonic  hemispherical 
pressure  distribution. 

It  was  next  decided  to  see  how  simply  this  distribution  could  be 
approximated.  Of  the  numerous  methods  for  approximating  the  pressure 
distribution  (power  series,  trigonometric  series,  etc.),  it  was  felt  that  a 
simple  relation  of  the  form 


P  =  Ps  cos"  9  (10) 

would  best  fit  the  purposes  of  the  study.  This  is  a  straightforward  choice 
since  this  relation  is  ox  the  familiar  Newtonian  approximation  form  and  requires 
only  the  generalization  that  the  cosine  exponent,  n,  is  not  assumed  to  be  2.  0. 

In  Figure  3,  P/pa  averaged  from  the  numerical  solutions  is  plotted 
versus  cos  S  on  log-luy  paper.  Thus,  it  P/pg  =  coo“  6  ,  we  can  write 
log  P/ps  =  n  cos  ®  and  n  he  the  slope  of  the  curve.  Fairing  a  straight 

line  through  these  data  produces  a  resulting  slope  of  n  =  2.  24.  Lines  which 
bound  the  whole  curve  are  seen  to  fall  between  2.  14  and  2.  34.  Thus,  a 
maximum  error  in  n  of  z -4  is  assumed  to  exist.  The  curve 

F/Ps  =  cos^’^40  is  also  shown  plotted  in  Figure  Z  for  comparison  to  the 
numerical  and  experimental  results. 

c.  Attitude  Equations  at  Hypersonic  Mach  Numbers 

The  pressure  distribution  selected  for  the  hemisphere  probe  was  seen 
to  be  of  the  form 


P  *  Pj  COS  8  ,  r\  •  2.Z4  (10) 

where  &  is  the  angle  between  the  stagnation  point  and  some  orifice  at  which 
pressure  P  is  measured.  With  reference  to  the  body  axis  system  defined  in 
Section  II.  1.  ,  the  attitude  angles  ©C  and  ^3  can  hi  determined  from  this 
assumed  pressure  distribution  and  from  the  pressure  measurements  in  the 
x-z  and  x-y  planes,  respectively. 
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Since  the  location  of  any  orifice  can  be  measured  in  terms  of  its 
angular  displacement  irom  the  stagnation  point,  herein  called  the  angle  6  , 
(or  any  orifice,  i,  the  pressure,  Pj  can  be  expressed  an  a  function  of  the 
stagnation  pressure,  P8,  and  the  cosine  of  the  angle  6  . 

Pi  -  Ps  cos'1  0  (1 

For  an  orifice  in  the  x-z  plane,  (see  Figure  1)  say  Pj 


p,  - 

Ps  cos" 

cos  8  * 

u/  U« 

COS  CL  “ 

a/(u^ • 

cos  8  m 

K-h  ur* 

cos  0  - 

cos  CL  c 

p,  = 

n  h 

Ps  cos  < 

ful  relation, 

since  ^ 

so  that 


in  the  x- z  plane;  thus  the  pressure  at  any  attitude  location  in  this  plane  is  independent 
of  <5  .  Any  pressure  in  the  x-z  plane  is  thus 

Pi  ■*  Ps  co»h  ci  cos  (oL  -  }  (17) 

where  is  the  orifice  angular  displacement  from  the  x-axis,  in  the  x-z  plane. 

The  determination  of  sideslip  angle  »  as  defined,  is  identically 
the  same  problem  as  the  determination  of  OL  .  All  that  is  required  to  show 
this  is  the  introduction  and  subsequent  cancellation  of  an  additional  angle,  say 

e  :  ,  Jfc 

.  -I  W 


s  i  n.  -=s=- 
U. 


such  that  the  pressure  distribution  in  the  x-y  plane  is 

Pj  •>  Ps  cos”  S  *  P5  cos  *  (j3  ■*  0j)  cos"  t  (19) 

As  can  be  seen  from  Eq.  (19),  the  distribution  in  the  x-y  plane  is  independent 
of  £  as  was  the  pressure  distribution  in  the  x-z  plane  independent  of  c> 

(1)  Determination  ox  Angle  of  Attack 

The  pressure  on  the  hemisphere  at  some  point  in  the  x-z  plane 
was  seen  to  be 

P  P*  cos"  ©  »  P»  cos  cos"  (ol  -  4>)  (17) 
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where  tlic  angle  &  is  invariant  in  this  plane.  Thus,  tin  pressure  distribution 


in  the  x  z  plane  is  independent  of  P  ami  the  angle  o.  Similarly,  an  expression 


for  CL  can  be  nondimcnsional-zcd  to  make  it  independent  of  Ps  cosn*f. 

Consider  the  pressure  a!  three  orifice  locations  in  the  x-z  plane  (sec  sketch,  p.  q). 


p,  -  p„ 

cos”  & 

*  / 

cos  (ol  - 

0,) 

(17) 

pt  »  p. 

co^S 

n 

COS  OL 

,  =  O 

(20) 

*» 

CL 

1 

o? 

cos  nS 

cos  (a.  - 

■0s) 

(21) 

• 

ro  all  three 

equations.  Thus, 

Equations  1  7, 

20  and  21 

represent  only  two  independent  equations  for  the  determination  of  unknowns 
(Ps  cosn<5)and  (X.  ,  or  three  independent  equations  for  the  determination  of 

(1  8  cosn<$),  OC  ,  and  n  if  it  is  not  assumed  that  n  is  known.  For  the  hypersonic 
ci  se  where  n  is  known,  any  two  of  the  three  eauations  are  sufficient  for  tiie 
solutions  of  these  two  variables.  However,  the  solution  for  OL  can  make  use 
of  all  three  equations,  if  the  inclusion  of  the  third  equation,  although  redundant, 
will  decrease:  the  uncertainty  of  OL  for  given  uncertainties  in  the  individual 
pressure  measurements.  In  general,  then,  the  choice  of  two  or  tin  ec  pressure 
measurements  will  depend  on  accuracy  rather  than  need. 


a.  Angle  of  Attack  Determined  from  Two  Piessurcs 


Equations  l7  and  20  can  be  used  to  eliminate  cos11  5  by 
either  cf  the  two  formulations: 


P.  -  Pt,  _  cos  YtX-0,)  -  COS* CL 
P,  -+  Pz  Co (oi  -  0 +  coshcC 


_ R _  cos  n  (<x  -  0>) 

P,  cos’’1  CL  (2  3) 

Equation  (22)  is  shown  plotted  in  Figure  4  for  values  of  0, 

chosen  to  cover  the  angle  of  attack  ranee,  remembering  from  Figure  3  that 

•  D  _  n  ”  n  .  i-7  .  i  c  _  ..  ,  .  i  i 


the  pressure  distribution,  P  ■  P*  cos  ©  ,  is  valid  only  for  ©  up  to  about 
65*. 


/  p_p  >  The  curves  of  Figure  4  demonstrate  the  sensitivity  of  the 

ratio(  p:  +  jto  CL  lor  different  values  of  (fi)  and  can  thus  be  used  to 

select  0,  's  that  will  minimize  uncertainties  in  OL  due  to  uncertainties  in 
measurements  of  Pj  and  Pj.  First  of  all,  however,  it  is  necessary  to  find 
an  express^  n  for  the  change  in  ( ~  due  to  changes  in  P,  and  p£.  This 

is  done  by  differentiating  Equation  (22)  ; 


i  \-£7TPTj 


'  R-P* 

Th 


To  be  complete,  it  should  be  considered  that  n  aiso  can  have  an  uncertainty 
T)~  A  n  ;  this  ui. certainty  will  be  considered  in  the  final  error  uncertainty 
equations.  " 


Performing  the  differentiation  in  Equation  (24)  and  dividing  through  by  Equation 
(22),  results  in 


z  Px  <LPX 


Z  P,  dPz 

(p.-p it)(pi^7 


(25) 


etc. 


.  A R  AF»,  +  s~>  ,  .  , 

Let  ~:j  ~  — =p“-  c  and  consider 

'  i  .  •*_ 


Equation  (25)  then  reduces  to 


■4C  P,  Pi _ 

(p,  -  P;.)  (p.  *  Pi) 


(26) 


For  various  values  of  OL  and  for  the  case  of  one  percent  pressure  errors, 
C  -  0. 01,  values  of 


A 


pL-.EkJ\ 
V  p.  -  P*  ) 


4C  P,  Px 

(p,  + 


(27) 


were  calculated.  These  increments  were  then  added  and  subtracted  to 
in  Figure  4  to  obtain  values  of  £  A  OC  (see  sketch  below). 


These  values  of  iAoC,  obtained  giaphicallv,  wei  e  then 
plotted  versus  OC  for  the  different  0  1 9  in  Figure  5.  As  can  be  seen,  the 
best  choice  of  0,  that  will  cover  the  OC  range  fre^m  -20"  to  50°  and  maintain 
(ac-<p)  iC  65°  is  0,  =  45*.  Foi  thia  case,  for  1%  errors  in  pressure 

over  the  entire  &-  range. 

Next,  consider  the  alternate  formulation  of  Equations  (17)  and 
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**m&mj*  fc  j  i*  i » 


(20), 


'fr  '  > 


l_  -  -  ‘  l  ,~.v  w:- 

■■•'  r: 


R 


ri  /  /  \ 

CQ5  (oC-  0.,/ 
C^S  n  OL 


(23) 


Equation  (23)  is  platted  in  Figure  6  fci  (he  0,  values  considered  in  Figure  4. 
The  uncertainty  in  OC  due  to  uncertainties  in  Pj  and  was  determined  in 
the  same  manner  as  was  done  for  Equation  (22). 


and  as  before. 


tfeA) 


.  idr 


ctR. 

P, 


dh. 


(28) 


<A- 


iR 


R 

A  (R/Pz) 


or 


(-R 


Pz  ; 


<tP2 

P2 

AP, 

P, 

2C 


aP2 


=  ±  C 

=  2.C. 


A 

P, 


(29) 

(30) 


Using  Equation  (30),  graphical  values  of  A  CL  were  obtained 
from  Figure  6  and  are  shown  plotted  in  Figure  7.  Comparison  of  Figures  5 
and  7  shows  the  results  to  be  practically  equivalent.  Thus,  Equation  (22) 
presents  no  advantage  in  accuracy  ever  Equation  (23)  and  the  simple  relation 


P, 


cos 


(oL-  0,  ) 


Pz 


\  -  -»/ 


appears  to  be  the  better  choice  for  angle  of  attack  determination  for  the  two 
pressure  measurements.  Furtliermoi  e,  (Z>.  -  45°  is  again  the  best  choice  of 
orifice  locations.  For  this  value  of  0(  ,  Equation  (.23)  reduces  to 


R 


p,  *  (tt)  0  ♦ta"°0  • 


which  can  nov/  be  solved  d.rectiy  for  angle  of  attack. 


CL  -  tan 


-l 


- ' 


(31) 


(32) 


Equation  (32),  therefore,  repiesents  a  closed  form  solution 
for  cX  using  o^ly  two  pressure  inputs  in  the  x-z  plane.  Since  this  equation 
results  in  r.o  loss  in  accuracy  over  the  alternate  orifice  calibration.  Equation 
(23),  and  has  the  advantage  of  a  closed  form  solution  in  CL  ,  it  has  been  selected  as 
the  most  promising  form  for  the  two-pressure  input  case. 
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■  T  ■  ■  ijr* 


An  analytical  expicssion  for  angle  ol  attack  uncertainty  can 
now  be  found  by  tak;ng  the  total  derivative  of  Equation  where  we  will  now 

allow  n  to  be  uncertain.  We  will  allow  0,  to  be  variable,  3ince  we  have  not 
yet  shown  that  0,  -  45°  is  an  optimum  location  when  n  is  allowed  to  vary. 
Equation  (23)  can  be  written  as 


v»  /' . , 

coz-  lC*-- 


-  (c05  +  iin  Ian  c-~) 


n  4 

f  ■  u?  i 

f  i 


The  total  derivative  of  both  sides  of  Equation  (33)  is  taken,  resulting  in 

n-l  ? 

ZCh-p-j  -*  n  (cos<fit  +  sin  Q \  sec.  at  cLx.  ^ 

-h  (cos  0,  +  Sin  0,t»v\oc)  ioge  (coS0,4  Sin0,  tancl)<in 

Solving  for  cLqL  produces 

r 

ZC  -  a  (cos  0,  +  si  n  0  fact)  dn  (35) 

Next,  assuming  accumulative  errors,  and  ctoL-  £  AoC  ,  etc. 

r 

£oC  «  (oc£0,+tavwrt-)“-^-  ZC  ^ loge  (cos 0,+ rhrxfttanoL)  An  (36) 


isquaticis  (3t>)  is  presented  in 


Figure  B  for  (p,  -  15 


- 1 -  ,  -  ,  -  it  -  «=> - *  r  *  - '  ~  » 

45"  for  C  -•  ^WP  of  *.  01  and  for  a  constant  value  of  An  of  ±.  ?.  for  the  pur¬ 
pose  of  selecting  orifice  location  0,  .  Again,  it  was  assumed  that  all  errors 
would  be  accumulative.  An  orifice  location  of  0,  -  45",  as  expected,  is  seen 
to  produce  the  lowest  uncertainty  in  A (X  .  Actually,  this  error  may  simply 
decrease  as  0.  increases,  however,  45*"  is  just  about  ar.  upper  limit  on  0,  if 
we  want  to  measure  to  minus  15°  or  20°  angle  of  attack.  Next,  the  error 
dependence  due  to  An  is  shown  in  Figure  9  for  0,  -  45*,  where  actual  values 
of  An  were  obtained  as  a  function  of  cX.  as  follows.  On  a  plot  of  the  actual 
hemispherical  pressure  distribution  (Figure  3),  lines  of  COS  ©  for  various 
n  's  were  drawn  to  intercept  the  curve  at  S  'a  up  to  60°.  Then  for  each  0  , 
values  of  An*=  /?.24~n  were  obtained.  These  values  were  then  used  in  Equation 
(36)  for  C  =  .  01  to  obtain  the  results  shown  in  Figure  9.  From  Figure  9,  it 
is  seer,  that,  up  to  about  45°,  the  A.OL  error  due  to  the  AP  and  An  uncertainties 
is  only  about  6/10  of  a  degree,  which  is  felt  to  be  quite:  good.  Furthermore, 
ever  this  range  the  uncertainty  in  An  alone.  \a  negligible  as  can  be  seen  by  the 
dotted  curve  which  was  plotted  for  An  =0,  and  which  follows  the  solid  curve, 
up  to  45°.  Finally,  for  comparison,  the  graphically  determined  values  of  Acx. 
from  Figure  7  are  also  shown  replotted  on  Figure  9;  and  it  is  seen  that  good 
agreement  of  the  two  methods  of  Aot.  determination  is  achieved. 


i  1 
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b.  Angle  of  Attack  Determined  from  Three  Pressures 

In  using  Equations  (17),  (/.0)  and  (21)  to  determine  OC  ,  it 
musf  be  remembered  that  the  assumed  pressure  distribution  is  valid  only  for 
about  6  —  65°,  or  (C<-  0)  values  somewhat  less  than  65*.  Since  CL  can  be 

as  large  as  50°,  0  can  be  no  less  than  about  -15":  therefore,  thetotal  orifice 
spacing  can  be  practically  no  greater  thu.-i  the  spacing  chosen  for  the  two 
orifice  case. 

The  three  orifice  calibration  of  oL  ,  however,  can  bo  co  npared 
to  the  two  orifice  calibration  on  the  basis  of  accuracy  arc:  the  method  of  pressure 
measurement.  The  formulation  chosen 


Pt  -  P?  cos"  -  COs”(<cl-0  *) 

Pt"P3  "  co5"  (<X-<pz)  -  COS*  (at  - 


(37) 


incorporates  the  use  of  prof  sure  differentials,  as  opposed  to  absolute  pressure 
measurements,  which  may  be  an  advantage-  from  an  instrumentation  standpoint. 
To  cover  the  OL  range,  the  orifice  locations  were  selected  tc  be  0,  -  45°, 

0i  =  0*.  03  -15*.  For  03  =  -15*,  the  pressure  distribution  is  valid  up  to 

OL.  =  50 *,  however,  a  singularity  (i.  e.  ,  — — ►O)  in  the  function  as 

CL. — —  74. limits  the  useful  negative  range.  Equation  (37)  is  shown  plotted 
in  Figure  10  where,  for  the  0's  selected,  the  function  reduces  tc 


R  -  P7. _ I  { I  +  U«  &)_  __ 

Pi-  P3  "  I  -  .7,51 of)7'- 

The  error  equation  for  this  relation  turns  out  to  be 


(36) 


rj/  --  f" 


~k  ' 

-t: 


BIT 

Pt  ' 


er 


-i 


(39) 


Pa/ 


c  rjnA-  + 

1 1  p.-p, 


P«  ±  P>  1  P-Ps 


(40) 


The  uncertainty  in  OL  ,  due  to  accumulative  C  ~  ±  0.  Cl 
uncertainties  in  individual  pressure  measurements,  is  shown  in  Figure  1  i. 

As  can  be  seen,  the  ACX.  uncertainties  are  larger  than  the  two  orifice 
formulation  for  most  of  the  angle  of  attack  range  (compare  Figures  7  and  1  1). 
From  the  standpoint  of  accuracy,  therefore,  the  three  orifice  relation  appears 
to  be  les  r  desirable  than  the  two  orifice  relation. 
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(2)  Determination  o£  Angle  of  Sideslip 


Thu3  far  the  attitude  angle  equations,  orifice  selection,  error 
equations,  etc.  hat  e  been  evaluated  for  angle  of  attack  determination.  As 
pointed  out  in  Section  II.  1.  ,  however,  there  is  no  fundamental  difference 
between  the  determination  of  angle  of  attack,  OC  ,or  angle  of  sideslip,  . 

Thus,  all  discussion  thus  far  written  in  Section  U  2.c.(l)  forCLorifice  selection, 
etc.  i3  equally  apj  licable  for  x3 


For  the  angle  of  adark  equation,  a  (£> x  orifice  located  at  45°  was 
seen  to  give  the  beat  R  / relation  for  small  <pc  values.  Therefore,  although 
the  range  is  smaller  than  the  angle  oi  attack  range,  a  location  of  (f ,  =45° 

will  similarly  be  the  best  choice  for  the  ^3  equation.  Furthermore,  the  same 
arguments  that  showed  the  simple  pressure  ratio  relation  to  be  optimum  holds 
equally  well  for  /3  ■  The  equation  is  then,  for  ^  ,  =  45*  in  the  x-y  plane 


z3 


and  the  error  equation  is 


"tan 


-v 


(41) 


* -  (cot  ar,  +  t»ys)  to  [oge(co*i  +sm  2f(ta»yS)J  An 


(42) 


d.  Attitude  Equations  at  Large  Angles  of  Attack 

The  attitude  angle  equations  were  derived  for  angle  of  attack  usage  up 
to  approximately  50".  Thus,  the  choice  of  exponent  n  in 


COS*  0 


(10) 


was  selected  to  minimize  errors  in  CL  over  this  range.  It  is  o :  some  interest 
to  investigate  the  errors  in  OC  at  angles  as  large  as  85®.  This  was  done  simply 
by  extending  Figure  3  to  large  0  's,  Figure  12,  and  by  using  Equation  (36)  to 
generate  Figure  13.  Values  of  An  used  in  Equation  (36)  were  obtained  in  the 
same  manner  as  was  done  to  generate  Fi.  ure  9. 


From  Figure  13,  it  is  seen  that,  above  about  50®,  the  angle  of  attack 
error  AcC  increases  rapidly  to  rather  large  errors  ~6e.  However,  the  error 
in  OC  ,  although  large,  would  still  allow  use  of  the  attitude  equation 


OC  =  tan-'  \{Z 

\  ru  i  ■  i 

'0,  ■•45* 

for  angles  of  attack  as  large  as  85*  with  no  discontinuities  in  the  equation. 


Yn 


w  - 


(32) 


3.  ATTITUDE  ANGLE  DETERMINATION  OVER  THE  ENTIRE  MACH  I'UMrtK 
RANGE 

The  ba'ic  expression  used  in  the  derivation  ol  the  attitude  equations  at 
hypersonic  Mach  numbers  ir.  Section  H.  2.  c.  , 

p  =  ps  cos"  e  (io) 

was  found  to  held  with  sufficient  accuracy  with  respect  to  both  experimental 
data  and  theoretical  solutions  (Seciion  II.  2.  b.  ).  Furthermore,  it  was  found 
that  the  exponent  n  was  equal  to  approximately  2.  24.  If  Equation  (10J  is  also 
satisfied  at  lower  Mach  numbers,  nose  pressures  alone  could  continue  to  be 
used  to  obtain  simple  attitude  angle  expressions.  To  examine  the  validity  of 
this  pressure  distribution  relation  at  Mach  numbers  less  than  6,  data  compiled 
in  References  (15)  and  (2)  from  the  literature  are  shown  in  Figures  14  and  15. 

These  plots  present  %  versus  cosS  on  lo  i  scales;  thus,  a  linear  relation 
in  would  indicate  a  distribution  of  the  form ^'^•cas'O  where  the  slope  of  the 

line  would  be  equal  to  n.  For  the  most  part.  Figures  14  and  1  5  do  indeed  show 
that  a  cosine  relation  provides  a  satisfactory  approximation  to  the  pressure 
distribution  for  Mach  numbers  as  low  as  =  .  5.  The  slopes  of  these  curves 
(  H  ),  however,  are  not  constant  but  are  Mach  number  dependent.  Values  of  n 
measured  from  Figures  14  and  15  are  plotted  as  a  function  of  M^in  Figure  16. 
The  results  of  n  correlate  with  M®  quite  nicely.  Furthermore,  an  extrapolation 
of  the  curve  should  satisfy  the  value  of  n  (2.  24)  assumed  to  hold  at  hypersonic 
Mach  numbers,  which  it  appears  to  do  reasonably  well.  This  result  is  in 
good  agreement  with  Reference  (15)  which  showed  a  cosine  exponent  to  vaiy 
from  1.  5  to  2.  3  in  going  from  low  supersonic  to  hypersonic  Mach  numbers;  in 
that  report,  the  exponent  was  defined  through 

P  —  "5"  R  f  *  -I  \ 

- i —  “  CU*>  « \J  (l3) 

p _ Lp 

*  3t  ra> 


This  equation  differs  slightly  from  the  present  correlation  having  the  disad¬ 
vantage  of  requiring  a  knowledge  of  free  stream  static  pressure  for  attitude 
determination. 


If  Mach  number  was  a  Known  input.  Figure  16  could  be  used  to  obtain  n 
accurately  enough  to  use  the  attitude  expression, 

K. 

Oi  *  tan"1  [Vr(-K  )  -  ij  (32) 

Without  being  given  ,  however,  an  expression  for  is  needed  which  is 

dependent  on  input  nose  pressure  but  independent  oi  n.  A  relation  of  this  form 
will  now  he  developed. 


Consider  again  the  three  basic  pressure  equations  that  can  be  written  for 
three  pressure  inputs  in  the  x-z  (pitch)  plane: 

»  Reas"  6  =■  F*  C0>r  <D  cos"  (oL-  (f>,)  (17) 


n  o  r  H  , 

-  Ps  cos  o  cos  cl 


Rj  =  PscosnS  cos"  (<X-^3) 


Pressures  Pi  ,  P2  and  Pjh  are  measured  quantities  and  it  is  desired  to  find 
OL  .  The  quantities  Cos  &)  ,  and  h  are  also  unknown,  but  three  indept  ndent 
equations  exist.  First,  (P5  cos*'§)  is  eliminated  by  dividing  Equation  (17)  by 
Equation  (20)  M 

-id-  3  L^-fe^li.1  (23) 

Pz  ^  COS  cC  J 

and  Equation  (21)  by  Equation  (20). 

-Ea.  »  [-  (44) 

P*  I  COSOL 

Next,  In  can  be  eliminated  by  taking  the  log  of  Equations  (23)  and  (44), 


'“5  (i)  ■  n 

•°*  (i-)  - "  * 


and  then  div.ding  Equation  (45)  by  Equation  (46). 


lag  (g/% )  =  iorf  Cos-0 

,os  (co^-^/coscl) 


Equation  (47)  expresses  CL  as  a  function  of  only  P,  ,  rz  and  P3  . 

We  now  have  an  expression  for  angle  of  attack  which  is  independent  of  y\ 
and  requires  only  a  pressure  distribution  of  the  form  P=f^coxM0  ,  regardless 
of  the  value  of  n  .  A  plot  of  this  function  for  0.  ~  45p,  03  =  -15°  is  given 
in  Figure  17.  For  the  case  of  0,  =  45°,  03  =  -45*,  Equation  (47)  reduces  to 


loaQypJ  .  PT  (i  4-tan  oc) 

•°£  (%)  lo<g  ~Y~  0  ~ t&lr'  <*0 


(47a) 


Using  Equation  (47),  it  is  possible,  therefore,  to  obtain  ot  (or^  )  independently 
of  n  .  Once  OL  (and  W  )  are  found,  however,  it  is  possible  to  obtain  h  by 
using  Equation  (32).  Taking  the  log  ox  both  sides  of  Equation  (32)  and  solving 


using  Equation  (32) 
for  r\  one  obtains 


lotf 

loS  (  vTT  J 

The  necessity  of  knowing  h  even  when  OL  and  ft  have  already  been  determi 
will  be  seen  in  Section  Ill  wherein  the  calculation  of  stagnation  pressure,  F s 
is  required  in  order  to  obtain  true  and  indicated  air  speed. 
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mined 


Equation  (47)  would  appear  to  be  very  useful  where  n  varies  but 

is  not  k  iown  as  an  input.  The  usefulness  of  Equation  (47)  at  hypersonic  Mach 
numbers  should  alsu  be  considered.  Obviously,  an  equation  that  can  be  used 
continuously  over  the  whole  Mach  number  range  is  more  desirable  than  having 
separate  equations  which  would  cover  incremental  Mach  number  ranges.  A 
Caparison  of  lib  nation  Equations  (47)  and  (32)  should  consider  the  following: 

AP 

1)  Accuracy  in  <X  (or )  for  given  ~p~  errors. 

2)  Range  of  use  in  oC 

3)  Simplicity  of  attitude  expression. 

OLviously,  Equation  (32)  is  the  simpler,  but  to  make  a  valid  comparison 
of  items  1)  and  2),  it  is  necessary  to  first  attempt  to  optimize  Equation  (47)  for 
accuracy  and  range  of  use  for  various  orifice  locations,  as  was  done  for 
Equation  (37),  A  series  of  five  different  or.  ice  locations  were  considered 
and  are  shown  in  Figurej  17-21. 


Before  these  figures  arc  discussed,  it  is  useful  to  consider  first  the  errors 
in  angle  of  attack  which  could  result  from  errors  in  pressure  measurements 
using  Equation  :7).  Taking  the  total  derivative  of  the  left-hand  side  of 
Equation  (47) 


dm  . _ L _ Ur  .  as] .  (x)  fas  _  *5.1  « 

I. p'  P*J  [*  (-StT  1. p’  p>  J 

JLP  A  P 

for  — pr~  =  — pr4-  ,  etc.  Equation  (49)  can  also  be  written  in  terms  of  the 
angle  relations 


Aftp) _ i — ^(aE.asl- Jmu 

"MrSj£*2)L p  P*J 

<•  /  > 

Dividing  Equation  (49)  Jlr/  results  in 


AR  A§1  I.r  [as  ap, 


{ cojfoi-  <3j)V|  * 


(49b) 


Next,  for  the  right-hand  side  of  Equation  (47), 


(50a) 


Equating  Equations  (50  a)  and  (51),  i.  e.  , 


I  (tavicl-tanCot-^j)) 


A-f  CP)  _  Af  (oe) 


(51) 


and 
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solving  I01A0C  results  in 


Act 


lot}  (coS^jts'uiQ^txn  ■flog  (cos  jZ^+  siw^.tsuri)  |S*| _ 

h  Cco*  sin  ^jtsh^Ctsnof-tasC^-  ♦Vn^ t»»w)(Uk<c-tkh(<- 


(52) 


The  numerator  of  Equation  (52)  is  of  the  form 


(53) 


and  it  was  found  arithmetically  that,  when  A  and  B  have  opposite  signs,  an 
error  of  -"-A P*  /'A+B^  could  occur,  whereas  when  A  and  B  have  like  signs, 

P  V  )  2Ap  ZAP 

the  maximum  possible  error  is  — p —  A  or  ~p —  B  »  depending  on  whichever 
of  A  or  B  was  numerically  greater.  Shown  in  Figures  22-26  are  plots  of 


AcL  versus  OL  prepared  using  Equation  (52)  for  AP/f3  =  .01  fui  the  five 
orifice  locations  of  Figures  17-21.  It  was  found,  however,  that  when  Af(P) 
was  large,  Equation  (52)  gave  poor  estimates  of  ±A OL.  This  resulted  because 
the  equation  assumes  a  straight  line  variation  in  (P)  over  the  corresponding 
AO C  interval,  which  was  not  always  the  case.  Therefore,  for  large  &yp  , 
values  of  A-f(P)  from  Equation  (49a)  were  calculated  and  graphical  values  of 
AOL  were  obtained  directly  from  Figures  17-21.  These  values  are  also 
shown  in  Figures  22-26. 


The  optimization  of  Equation  (47)  can  now  be  obtained  through  examination 
of  Figures  17-26.  The  results  of  these  plots  are  given  in  Table  I.  which 
comprises  a  summary  of  the  investigation  listing  cases,  accuracy  uncertainty 
due  to  1%  pressure  uncertainties  at  C*-  =  0  and  range  ol  use.  The  purpose  of 
examining  the  five  different  ca3es  is,  of  course,  an  attempt  to  achieve  for 
Equation  (47)  a  calibration  which  will  produce  a  maximum  attitude  angle  range 
of  use  and  a  minimum  sensitivity  to  pressure  error.  These  requirements  are 
somewhat  conflicting,  however,  as  can  be  seen  in  Table  I. 


Table  I 

ANGLE-OF-ATTACK  CALIBRATIONS 
FOR  LOGARITHMIC  RELATION 


POSSIBLE 

LIMITS  IN  RANGE  OF  00 

0, 

<*3 

INACCURACY 

•  UPPER 

LONER 

GURE 

CASE 

ATOC-0* 

MIN. IN  CURV 

Z«r  0,-0  | 

17 

1 

45* 

-15* 

14*.  -3.0* 

62* 

60* 

-  7.5 

-20 

18 

2 

15* 

-30* 

i . 9* ,  -1.6* 

45* 

35* 

-15 

-20 

18 

3 

45* 

-45* 

0.5*,  -1,0* 

38* 

20* 

-22.6 

-20 

20 

8 

80* 

-30* 

2.1*.  -1.7* 

49.5* 

35* 

-15 

-  5 

21 

6 

30* 

-30* 

1.7*,  -1.4* 

42.5* 

35* 

-15 

-35 
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An  explanation  of  the  limits  in  attitude  &ngk 
given  as  follows: 


ranee  sriOwn  in  iudic  i  arc 


1)  Minimum  in  Curve  -  Causes  an  upper  limit  in  attitude  range 
because  an  attempt  to  use  the  calibration  above  this  minimum 
would  result  in  a  double  value  attitude  solution  for  a  given  set 
of  pressure  inputs.  A  cross  plot  of  the  maximum  positive 
range  in  QL  limited  by  a  minimum  in  the  calibration  curve  is 
shown  in  Figure  27. 

2)  |oL~0(  |  =  65  e  or  |cX— 03|  =  65“  -  The  calibration  equation  was 
uerived  from  an  assumed  pressure  distribution  of  P  =  fjC05  0 
which  is  considered  to  be  valid  only  up  to  an  angle  between  the 
orifice  and  the  stagnation  point  of  approximately  05°. 

3)  (ZoC—  )  =  0  -  The  functions  become  discontinuous  if  this 
occurs,  see  Equation  (47). 

The  range  in  crificc  location  was  selected  in  view  oi  these  limits  and 
desired  range  of  angle  of  attack  (~ZO°  —  OC.  ^  50°)  as  follows: 


1) 


2) 


3) 


Upper  limit  in  0,  was  set  by  condition  2  above  (for  cos  Q 
distribution  to  hold  to  -20°,  0(  should  be  no  greater  than  about 

45“;. 

Upper  limit  in  l&i  was  set  by'  condition  1  (upper  CL  limit  in 
calibration  curve  due  to  minimum  decreases  as  |03|  increases) 
or,  condition  2  (for  cost)  distribution  to  hold  to  50°,  W  should 
be  no  greater  than  15°). 


Lower  limit  in  |03j  was  6et  in  consideration  of  the  fact  that 
inaccuracy  in  (X.  near  CX.  -  0  increases  as  |0oi  decreases,  _ 
by  condition  3  (for  oc  range  to  -20°,  |0j|  should  be  no  less 
than  40*). 


the 

cUivi 


Considering  the  possible  OC  inaccuracies  at  CX  =  0  listed  in  Table  I,  it 
can  be  seen  that  must  have  a  large  negative  value,  perhaps  45°,  to  make 
the  calibration  of  any  practical  use.  This  requirement,  however,  restricts 
the  upper  limit  in  useful  OC  range  to  about  +20°.  Thus,  in  terms  of  accuracy 
and  range  of  use,  Equation  (47)  is  not  as  acceptable  as  Equation  (32)  in  the 
hypersonic  range. 

In  addition,  it  was  previously  felt  that,  even  though  the  required  OC.  range 
was  -20"  to  +50°,  it  would  be  most  desirable  to  have  the  calibration  operable 
at  large  limits  in  oc.  (see  Section  II.  2.  d.  ).  Such  is  not  the  case  for  Equation 
(47),  since  discontinuities  in  the  calibration  occur  at  negative  OC  values 
(within  the  required  range  for  all  but  03  =  -45°,  case  3)  and  multiple  values 
occur  in  the  positive  CX.  range. 

One  possible  solution  or  compromise  would  be  to  use  one  equation  as  a 
back  up  for  the  other.  At  high  altitude,  high  Mach  numbers  when  the  attitude 
angles  could  conceivably  be  quite  large.  Equation  (32)  would  be  most  valuable. 
A  switch  over  to  Equation  (4  7)  at  supersonic -transonic  Mach  numbers,  where 
the  attitude  angles  would  not  be  too  great  (say  —  20*)  could  then  be  made. 

19 


■V. 


Orifice  locations  0, ,-  45s.  ($.  =  0s  arc  romnatibk  for 

additional  orifice  input  could  be  made  available  as 


both  solutions  and  the 
needed. 


The  required  attitude  range  in  is  far  less  demanding  than  for  CX  ,  being 
from  -15s  to  +15°.  Thus,  there  is  no  requirement  to  make  tft  or  ,  the 
orifice  locations  in  the  x-y  plane  (see  sketch)  small.  The  J  ang'es  should  be 
greater  than  30s  tc  eliminate  any  discontinuity  in  the  calibration  and  could  perhaps 
reach  40°  or  45s;  they  are  limited,  however,  by  thecosO  distribution  at  large 
Oc  .  *  Actually,  o,  =  45  s,  &s  =  -45s  provide-  good  accuracy  ncai/3=  0,  do 
not  exceed  the  acceptable  cos  ©  range,  except  at  very  large  angle  of  attack, 
and  allow  the  simplified  calibration 


0 


"tan 


_l 


(41) 


for  the  hypersonic  case. 

Selection  of  the  logarithmic  equation  for ^ 

'■*  (-rr) 

in  place  of  Equation  (41  )  in  the  hypersonic  case  is  dependent  upon  the  same 
question  of  accuracy  and  convenience  as  was  just  discussed  lor  angle  of  attack. 
For  the  case  of  =  45s,  J 3  =  -45s,  Equation  (54)  reduces  to 


(54) 


(~S~)  _  lo£  0  + 

io4$r)  lo*  4*  ('  -  Un/3) 


i 

i 


t 

) 


At  OC  ^  45s,  -  los,  the  total  angle  between  the  stagnation  point  and  the 

orifice  at  - is  about  65s  for  K3  =  40  s. 
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SEGTIOiS  ill 

TRUE  AIR  SPEED  AND  INDICATED  AIR  SPEED 


1.  CALCULATION  PROCEDURE 

In  Section  II,  it  was  shown  that  pressure  distributions  in  the  x-z  and  x-y 
pta  nes  arc:  sufficient  inputs  to  determine  the  hemispherical  probe  attitude 
angles,  Oc  and  / 3  .  Using  ex  and  / 3  ,  the  centerline  orifice  pressure 
measurement  and  an  assumed  pressure  distribution,  the  probe  stagnation 
pressure  can  be  determined.  For  a  given  free  stream  density,  the  stagnation 
pressure  can  then  be  used  to  obtain  free  stream  velocity  (true  air  speed) 
through  a  stagnation  pressure,  dynamic  pressure  relation. 

2,  CALCULATION  OF  STAGNATION  PRESSURE 

The  pressure  distribution  on  the  hemispherical  probe  surface  is  assumed 
to  be  of  tlio  form 


P  =  Ps  COS  0  (10) 

where  0  is  the  angle  between  the  stagnation  point  and  some  orifice  at  which 
pressure  P  is  measured,  see  Section  II.  It  was  also  seen  in  Section  II  that 
if  the  orifice  is  located  on  the  probe  axis  of  symmetry 

CO50-  COS  8  COS  OC  (15) 

where  angles  6  and  OC  are  also  defined  in  Section  II.  1.  Substituting  this 
relationship  into  Equation  (10)  and  solving  for  P«  results  in 

Fj  •»  Rj ;  / COS  &  CoS  OC  (55) 

where  Fz  is  the  pressure  measured  at  orilice  2. 

Next,  consider  the  relations  between  angles  OC  ,  /3  and  <5  as  given  in 

Section  U.  1 . 


5  in  S  - 


(56a)  sinOL-  -  »  (57a)  sin/3  ■  *£ —. 


(58a) 


/  * 

(56b)  COSOL* TT  (57b)  cos/S.  ~~ ■  --,■( 58b) 

V  (U*tw  )  *  '  (ur  -+  Xr*)'* 


tan  5  =  ~~j~ — j-ir-  (56c)  tan  OC  »  — ^ -  (57c)  t »*/?  = 

(VX  +  {JJr  )'z  u- 


(58c) 


u. 


(59) 


By  inspection  Equations  (5bc),  (57b)  ami  (5be)  combine  to  give 

tan  S  *  cosoC  tan  ^ 


8  c  tan-'  [coscxtan/3 

L  /  . 


(59a) 


Combining  Equations  (55)  and  (59a)  results  in 


_ a _ 

COS  (t*n  ‘  (COIOL. 


(60) 


This  equation  then  presents  a  closed  form  solution  for  the  pilot  pressure 
in  terms  of  the  pressure  measured  on  the  centerline  orifice,  the  attitude  angles 
OC  and  /3  and  thi  cosine  exponent  Vi  .  For  Mach  numbers  greater  than 
6,  it  was  shown  in  Section  II  that  Yt  =  2.  24;  that  is.  Equation  (60)  can  be  written 


for  the  hypersonic  case. 


Pz 


£(cosol)  Co$(tav."  a-osectay9))J 


Z.7.4 


(60a) 


For  transonic 'Supersonic  Mach  numbers.  Equation  (60)  can  be  solved  only 
after  obtaining  n  through  the  method  described  in  Section  II.  3.  Equation  (60a) 
can  be  written  entirely  in  terms  of  measured  pressures  for  the  case  where  the 
attitude  angles  are  expressed  as  a  function  of  two  orifice  pressures  and  n  is 
known.  For  hypersonic  Mach  numbers  then: 


i 


Substitution  of  Equation  (32)  and  (41)  into  (60a)  results  in 


(41) 


Rt _ _ 

'  jcov jtw ' (vr ( ')]  (vT(^-i)J]j^(60b) 

J 

Because  of  the  complexity  of  this  expression,  several  simplifying  small  angle 
assumptions  were  considered.  First,  lor  the  assumption 

tTavi  <>  s  S  (6l) 


rj 


cos  *  *4  -  j)]  COr  tan 
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Eqiictio.ii  (60b)  reduces  to 


R»=  jf  os  -.)]  m  (l)^  -  i)  <62) 

which  offers  a  significant  reduction  in  complexity  over  (60b). 

The  error  in  fs  U3ing  (62)  is  shown  in  Figure  28  for  ’s  of  0°,  5®  and 
10°.  As  can  be  seen,  the  error  introduced  in  the  worst  case  is  only  .  07%. 

Thus,  Equation  (62)  introduces  negligible  error  in  the  pitot  pressure  calculation. 
If  is  possible  to  simplify  the  stagnation  pressure  relation.  Equation  (62),  further 
by  the  additional  assumption 

tan  ^3  =  0  (63) 

Therefore,  ~tsy\  &  -  O  and  Equation  (62)  reduces  to 

p 

'  *  (64) 

The  error  in  Ps  usin'*  Equation  (64)  is  shown  in  Figure  29,  again  for  /Q  's 
of  0®,  5®  and  10°.  Errors  for  this  case  can  be  as  large  as  3.  5%.  In  view  of 
these  errors,  Equation  (62)  would  appear  to  be  the  best  choice  of  the  expression 
for  ,  however,  the  value  of  the  trade  in  accuracy  for  simplicity  can  only  be 
ultimately  judged  in  the  case  of  a  specific  application. 

3.  RELATION  BETWEEN  STAGNATION  AND  DYNAMIC  PRESSURE 


a.  Imperfect  or  Real  Gas  Solution  for  Super  sonic -Hyper  sonic  Stagnation 

Streamline 

For  an  imperfect  gas,  *'  e  following  flow  conservation  equations  are 
valid  for  the  flow  in  a  streamtube  which  goes  through  a  normal  shock  and  comes 
to  rest  at  the  stagnation  point  of  a  blunt  body  (velocities  are  given  with  respect 
to  a  body  fixed  coordinate  system). 

First,  because  the  process  is  one -dimensional  through  the  shock,  the 
mass  flow  rate  entering  the  shock  surface  per  unit  area  must  exactly  equal  the 
mass  flow  rate  leaving  the  shock.  The  mass  continuity'  equation  for  this  process 
is  thi  refo'-e 


2.U. 


/V  Ur 


where  the  subscripts  are  noted  in  the  sketch  below 


FLOW 


/ 


"  5t*g  hm  ion  POIHT 


- j  r  •’“‘^ 

shock  surface  must  be  zero, 
momentum  equation  is 


shock  v.-a VC  is  SiaiiOiiai  y,  ljic  ..cl  lull  (_■  On  Inc 
The  force  balance  or  the  conservation  of 


+  /°cPto  -  Py  +  f> r  (66) 

Finally,  for  an  adiabatic  process,  the  energy  equation  between  the 
shocked  gas  at  station  r  and  the  stagnation  point  3  is 

jUr!  =  k  (67) 

Tliese  three  conservation  equations  form  the  basis  of  the  stagnation 
streamline  pressure  calculation.  Note  that  at  this-  point,  no  gas  restrictions 
(i.  e,  ,  perfect,  incompressible,  etc.)  have  been  assumed.  Rewriting  Equation 
(66)  and  substituting  for  Ur  using  Equation  (65)  results  in: 


K 

=  P.+/3.U/ 

1 

> 

$ 

£ 

N 

(68) 

Pr 

- 

0-#) 

(68al 

The  thermodynamic. 

P 

equations 

r  Z^RT 

imperfect  gas  law 

(69) 

h 

=  CpT 

definition  of  specific  heat 

P 

at  constant  pressure” 

(70) 

are  next  substituted 

!  in  Equation  (67) 

resulting  i.n 

(S§~)  - 

.  iy! 

,  L£*-)  _B_ 

(71) 

vz*4 

/7,  2  r 
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.  Z  R  4  />, 

Equation  (VI)  is  next  solved  for  P»  ,  the  stagnation  pressure: 


(7U) 


Equation  (70)  merely  constitutes  a  definition  of  Cp  ,  it  is  not  meant  to  imply 
that  this  value  is  insensitive  to  ’“"perfect  gas  properties  (i,  e.  ,  a  constant). 


! 
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U'i“6  E’Ui*tio"  1681  «»  «““*■>*<«  Pr  ■  »nd  defining  dy.rtu.uc  pte..u,«, 

f  "  ~/°“  U“  <7 


(71b) 


(73a) 


b'  Sta8nati°n^namiC  P—  delation  for 

to  dynamfcqpfes°sar7eia^i"  ^  ^er£ect  S^^esult  for  the  relation  of  stagnation 
/_Cp_®  P  essure*  Vo.  *  However,  analytical  relations  for 

v  2-R  ^ r  must  be  formulated  to  solve  the  equation  resulting  in  =>  T  i*«--  c 

.o^rrto^rt;r»  T  ~lr 

F*/cn  can  be  found  fortfe  fr  ,h*C  calculated  exist  in  tne  literature  so  that 
a^nn.pnon,.  oornparn  «h?  result,  to  the  c-„ct  sAion,.  “ - „ 

vai  iauui!  nt-Lwopn  pvinj-  -»**-!  <-4^ 1:  r  .  j  -  .  .  -  r*  ' — 

functions  This  i,.K - solutions  witn  simple  correlation 

*  --  —*•••  <**«»* 

First,  assume  that  the  fluid  behaves  as  a  perfect  gas  allowing 


P  yRT, 

2r  *  Z5  -  f 

(74) 

A  _  cP 

r 

R  Cp  -  Cy. 

a7- 1 

(75) 

Suostitute  Equations  (74)  and  (75)  into  (73a)  which  gives 

**/  -*7  2  + 

8  /  r  1  ?' 

-  a  fz-  f  y 

A  \  m/ 

(76) 
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which  reduces  to 


/)  .  >'  !  .  T 


(?6a) 


Next,  assume  the  gas  between  stations  cc>  and  \f  is  incompressible, 


in  which  case 


f-  'f- 


Ll± 

1 


see  Relerence  (17) 


By  substituting  Equation  (77)  and  (78)  into  (76a) 


8/  =  2  +  -E=-  -  &- 

/{i  1  F 


For  hypersonic  conditions 


p*4  * 0 


Thus,  Equation  (79)  can  be  written 


Ps  /  =  2.  — 


Equation  (81)  becomes  recognizable  now  as  the  most  commonly  derived  solution 
for  the  stagnation  pressure. 

The  next  and  final  approximation  requires  taking  an  average  or  constant 
value  fcr  over  the  range  of  usage  and  thus  simply  assuming  a  constant  ratio 

of  .  r* 


n  /  . 

J  tr  /  -  «  K 

't 


The  results  of  the  solutions  are  now  shown  tor  the  velocity  envelope 
of  interest  in  Figure  30.  It  is  interesting  to  note,  that  the  assumptions 
and  incompressible  flow  have  had  little  effect  on  Equation  (76a),  thus,  0 

Equation  (81)  suffers  very  little  from  these  assumptions.  Logical  candidates 
for  an  approximation  to  Eq.  (73a)  thus  become  Eouation  (81)  or  (82). 


A  correction  to  Equation  (81)  or  (82)  could  include  velocity  and  alticude 
or  velocity  and  density,  e.  g. 

F  +  f(Fu")  (83’ 

U  ' 

However,  J  -  is  quite  insensitive  to  P &  and  the  correction  cat  be  written 

simply  in  terrm  of  U*,  ,  ' 

F*/7  =  1.84  +  -f(U»)  (83a) 


Finally,  front  a  c  irve  fit  of  Eq.  (83a)  to  Eq.  (73a)  for  the  range  of  interest 
(Fig.  30).  the  results  are  _  .  fl4  U«F  (83b) 

y<fr  ~  l,cs  4  moo 


26 


Equation  (8 3  b)  is  plotted  iii  Figure  30  and  is  seen  to  hoid  to  within  ±1/2%. 


4.  TRUE  AIR  SPEED,  HYPERSONIC 

Equation  (60n)  can  be  substituted  into  Equation  (83b)  to  express  true  ail 
spi  ed  in  terms  of  the  measured  quantities  OL  ,  jQ  ,  and  density.  This 
velocity  relation  would  still  need  to  be  salved  fbr  Ua,  ,  however,  which  can 
not  be  done  directly.  An  iterative  procedure  whereby  values  of  Uo>  would  be 
assumed  is  needed  to  obtain  a  solution  in  u*,.  Rather  than  resort  to  this 
procedure,  it  was  decided  to  first  investigate  the  necessity  of  the  complexity 
of  this  expres  .ion  for  ■  s/cl.  This  was  done  by  investigating  the  error 
introduced  by  the  use  cf  Equation  (82), 

Ps/^=  K  (82) 

where  K  =  1.  92  ±.  03  since  1.  89  <  R/4,  <  i.<?<;  over  the  range  of  interest,  see 
Figure  30.  Then  from  Equation  (72),vwe  can  write 

(J2  ,  — Lfk —  (8  4) 

K 

Differentiating  this  expression  results  in 

-2  LL  -  -  ~T  +  -TT  ^  -  “£r-  ^  /°°°  C8  5> 

f  op  K  f-'ao  y/o»  / 

Next,  divide  Equation  (85)  by  Equation  (84),  assume  ~  ^‘X/'\j 

etc.  and  again  assume  all  errors  are  accumulative;  the  result  is 

A  Uct,  _  _|_  f  A  K  +  //•«»  1  tCb\ 


From  Equation  (82) 


A  K,  _  .03 


A  1WIU  4— »  v±  to.  L  ±  V-/ 1 1  (Ufc;  tv 

within  5%,  -  ±.05. 


n  Po. 

<  j 

=  ±.015.  Assume  density  is  known  to 


In  order  to  get  s/fs  >  we  rnust  differentiate  Equation  (60),  which  woulc 
be  difficult  and  unwieldy,  but  it  can  be  done  for  the  simple  case  o>  ji  -■  0, 
Equation  (64),  which  will  still  serve  the  purpose  of  illustrating  the  validity  of 
the  s  K  assumption.  Equation  (64)  In  terms  of  OL  is  simply: 

R  “  '  *■/  cos*  oL 

and 

diP^  a  R  n  cosci  sin  of  dot  +  coi,nCLcl  Pz  +  R  coa*of  log  cos olJL>\  (88) 

Again  assume  d.  Ps  -  A  R  ,  etc.  giving  * 

R  b 1  n  Of.  /,  A  R _ ,  i„ ...  rr.urt  a  •«  /qq\ 


Cos  U. 


+- 1  o q  cos  of  tin 

R 


Substituting  this  expression  hack  into  Equation  (86),  the  result  obtained  is: 

^Lk=„L[4li  +  +  +  cosctAn")  (I 

U*,  2.  |^K  p*  cos’'7oc  Ht  j 


Note  that  from  Section  II,  for  (p ^  =  45° 


(36a) 


Therefore, 


If  this  expression  is  evaluated  at  OL  =  0  for  illustration  purposes 

1  1  2.  [_  K  Pa>  J 


(86c) 


AUlO 

U« 

AU» 

Ua> 


=  1/2  (.  015  +  .  05  +  .  01)  =  .  0375 
=  3-3/4% 


From  this  result,  we  can  see  that  the  error  in  U*  due  to  the  uncertainty- 
in  K  contributes  only  3/4%  to  the  total  3-3/4%  uncertainty.  Therefore,  the 
error  in  U,*,  introduced  by  the  assumption  of  the  simple  relationship  's/fe  ~  K 
is  small  and  we  are  justified  in  using  'yo.  =  1.  92  for  the  velocity,  density, 
stagnation  pressure  relation,  D 

Solving  Equation  (84)  for  gives 

u-'l02l( jk?  ,9' 


Substitution  of  Equation  (60b)  into  Equation  (90)  produces 

/  R  & 

1  !  1-021  (-^k) 

“  =  |cm[W(Vz(^)*-i)]  cos[W[c.,|U-'^)*-lI|(V2f^.)*-|J]]p  (91) 

which  gives  true  air  speed  in  terms  of  F?  ,  P*  ,  Rj  and  . 

Since,  however,  it  was  shown  in  Section  III.  2  that  the  assumption 

tan  5  =  5  (61) 

introduced  negligible  error  in  F$  ,  Equation  (91)  can  be  rewritten  for  M  =  2.  24, 
including  this  approximation,  as  i_ 
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Figure  31.  For  up  to  10°,  Equation  (92)  gives  very  good  results  over  the 
entire  (X.  range;  and,  therefore,  represents  a  satisfactory  approximation  to 
the  air-data  velocity  equation. 

Finally,  for  the  sake  of  completeness,  the  assumption  of 

tav>  {2>  =  O  (63 
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can  be  introduced  .nto  Equation  (92)  resulting  in  the  simplification 

,  ,  I OZI  (P*/pJ)^  , 

u-  •  ~  rt«--1/nVi»  -  Tmr~ 

The  error  introduced  in  U®  by  this  assumption  is  shown  in  Figure  32. 

5.  INDICATED  AIR  SPEED,  HYPERSONIC 

Indicated  air  speed  (j8  is  the  air  speed  obtained  using  sea  level  ambient 
density  reference  conditions.  For  the  methods  used  in  this  study,  Uj.  can  be 
expressed  simply  as 

M.  -  Lb(^)''4  <?-» 

whcreA  =  .  00237?  siugs/ft’l  The  value  of  U_f  can  be  obtained  once  U.*is 
known,  or  by  direct  substitution  of  Ecuatiop  (92)  into  Equation  (99)  which  gives 

i  i  _  1.02.1  /qc\ 


U  =  - 


cos 


Equation  (95)  avoids  introducing  the  uncertainty  of/Dp.  into  the  Uf  expression. 
The  simplifying  assumption,  Equation  (6l),  that  was  made  for  U®  is,  of 
course,  also  applicable  to  the  indicated  air  speed  expression,  Equation  (95). 


6.  AIR  SPEED,  SUPERSONIC 

Determination  of  the  air  data  outputs  U®  and  for  the  hemisphere  press  .re 
probe  is  dependent  upon  determination  of  stagnation  pressure  and  upon  the 
relationship  between  the  stagnation  pressure  behind  the  normal  shocK  is  and 
the  free  stream  dynamic  pressure  <2.  .  For  Mach  numbers  less  than  6,  it  is 
possible  to  determine  F$  once  n  isnound  as  discussed  in  Section  II. 

Fur  M®  >  6,  it  was  found  that 

K  (?2) 

TK 

where  K  =  1.92  i.  03  ever  the  hypersonic  range  (see  Section  HI. 4.  ).  The 
assumption  of  a  value  of  K  of  1.92  resulted  in  maximum  errors  in  U®  of  3/4%. 

For  6,  values  of  %  were  computed  for  Mach  numbers  as  low  as 

Mco  -  1,  anclthe  results  are  shown  in  Figure  33.  At  Mach  numbers  up  to 
about  M*  =  4,  's//<jp  decreases  with  increasing  M®  .  This  decrease  is  t.  ue  / 
Mach  number  dependence  existing  even  for  an  ideal  gas.  The  variation  in^/j^ 
in  the  hypersonic  regime,  bowevei,  (i.  e.  the  slight  increase  in  P«/ jP  with  M*  ) 
is  a  real  ga  s  effect  for  the  free  stream  velocities  and  densities  of  t.ne  re-entry 
trajectories  considered  in  this  study.  Errors  in  U®  resulting  from  the 
assumption 

=  1.92  (82) 
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are  shown  as  a  function  of  free  stream  Mach  number  in  Figure  3*1.  For  Mach 
numbers  as  low  as  M»  =  2,  the  true  air  speed  errors  resulting  from  the  use 
of  Equation  (82)  are  2%  or  less.  Therefore,  in  view  of  velocity  uncertainties 
of  ±2-1/2%  and  ±1/2%  resulting  from  assumed  uncertainties  of  ±5%  in  density 
and  ±i%  in  pressure,  respectively.  Equation  (82)  t.ppears  to  have  practical 
use  in  the  supersonic  range. 

For  Mach  numbers  less  than  ,2.  Equation  (82)  rapidly  becomes 

unusable:  for  example,  at  M*  =  1,  —  ~  20%.  Thus  it  appears  that 

the  hemisphere  pressure  inputs,  along  with  free  stream  density  and  speed  of 
sound,  are  inadequate  for  accurate  dote  rtnir.ation  of  true  air  speed,  indicated 
air  spied  and  Mach  number  at  M*.less  than  2.  Considerable  effort,  hov/ever, 
can  be  found  in  the  literature  related  to  the  use  of  a  hemisphere  cylinder  pitot- 
static  air  data  probe  at  subsonic,  transonic  and  low  supersonic  Mach  number; 
for  example,  References  (18)  and  (19).  Therefore,  based  on  these  references 
and  the  above  statement,  it  is  felt  that  the  pitot  static  tube,  that  i6,  the 
additional  input  of  static  pressure,  is  needed  in  order  to  determine  Mach 
number  for  2. 

With  density  given,  true  air  speed  (free  stream  velocity)  can  be  obtained 
for  M®  >  2  using  Equation  (82).  The  indicated  air  speed  calculation  merely 

requires  substitution  of  the  sea  level  density  for  A C ^  in  Equation  (82)  (see 

Section  III.  5.  ).  Therefore  '  / 

ux  *  2 
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SECTION  IV 


DETERMINATION  OF  ALTITUDE  AND  MACH  NUMBER 


1.  ALTITUDE  DETERMINATION 

An  altitude,  or  more  properly  a  pressure  altitude,  cannot  be  obtained 
from  the  pressures  measured  on  a  hemisphere  in  high  supersonic-hypersonic 
flow.  This  conclusion  is  a  direct  result  of  the  insensitivity  of  the  local  flow 
field  pressure  distribution  to  free  stream  Mach  number  and  hence  the  inability 
to  determine  M„  and  subsequently  static  pressure  from  the  input  of  surface 
pressure  (see  Reference  1  4  for  a  further  explanation  of  this  conclusion).  Thus, 
the  only  altitude  output  one  can  get  from  the  air  data  inputs  must  come  from  the 
input  of  density,  i.  e.  simply  a  density  altitude  (when  given).  Using  the  infor¬ 
mation  on  density  altitude  variations  available  in  the  1962  U.  S.  Standard 
Atmosphere,  Reference  (1^),  it  is  possible  to  estimate  how  accurately  one  can 
determine  altitude  when  density  is  known. 

A  detailed  investigation,  or  detailed  results  of  investigation.-,  of  atmospheric 
variations  and  uncertainties  in  pressure,  density,  etc.  will  not  be  attempted 
in  this  report.  This  information  is  adequ  ately  covered  in  Reference  (1  9)  and 
references  cited  in  that  work.  Rather,  we  will  look  at  the  conclusions  of 
Reference  (19)  only  to  the  point  where  we  may  formulate  the  ability  to  obtain 
altitude  when  density  is  given  as  an  air  data  input. 

Figure  35,  taken  from  Reference  (19),  presents  the  U.  S.  Standard  Density 
Altitude.  The  dotted  lines  in  this  figure  are  failings  through  density  extremes 
(symbols)  observed  at  given  altitudes.  From  these  extremes,  altitude  uncer¬ 
tainties  at  a  given  density  can  be  found  as  lying  horizontally  between  the  density 
uncertainty"  band.  Since  these  extremes  in  density  are  on  the  order  of  ±50%,  an 
additional  measured  density  uncertainty  of  ±5%,  which  is  to  be  assumed,  becomes 
of  little  practical  significance  in  determining  altitude  from  density.  One  finds 
that  the  density  is  very  sensitive  to  altitude.  Conversely,  altitude  is  a  weak 
function  of  density;  thus,  the  lar^e  mec-i  trtiulies  in  density  produce  much  smaller 
uncertainties  in  h  .  Values  of  *yV.  in  percent  are  given  in  Figure  36  as  a 
function  of  altitude. 

One  can  see  from  this  analysis  that  an  altitude  based  on  density  gives 
nominally  about  ±8%  uncertainty  in  h  .  Furthermore,  the  extreme  density 
une'ertainty  that  was  assumec  iO  exist  is  all  inc  usive  for  latitudinal  and 
seasonal  variations.  For  an  actual  flight  case  with  latitude  and  season  specified. 
Reference  (1  9)  points  out  that  the  uncertainty  in  density"  would  be  much  smaller. 
Using  appropriate  density-altitude  charts,  therefore,  altitude  could  be  determined 
to  much  better  than  ±8%  from  a  density  input. 

2.  MACH  NUMBER  DETERMINATION 

As  stated  previously,  the  pressure  distribution  on  the  hemisphere  and  the 
pitot  pressure  are  insufficient  inputs  to  obtain  in  the  hypersonic  regime. 

When  density  is  given,  however,  velocity  (true  air  speed)  is  calculable,  and 
Mach  number  can  then  ba  obtained  through 

(97) 

»•» 
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The  free  stream  speed  of  sound,  however,  must  be  obtained  from  the  density- 
altitude.  Thus,  we  must  look  at  the  uncertainty  in  &.«.  as  a  function  of  altitude. 
From  Figure  37,  also  taken  from  Reference  (19),  it  can  be  seen  that  at  a  given 
altitude  the  uncertainty  in  temperature  (symbols)  is  on  the  order  of  ±10%. 

The  variation  of  temperature  with  altitude  is  not  strong  enough  to  add  much 
additional  uncertainty  in  temperature.  Thus,  the  total  uncertainty  in  "T^  is 
still  on  the  order  of  ±10%.  Since 
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a  ±10%  uncertainty  in  X*  results  in  a  ±5%  uncertainty  in?t„. 


Now,  since 


AM„  _ 


U a, 
AUo 

UL 


A&a 


Mm 

In  Section  III,  it  was  shown  that  — ~  ±4%  so  that  the  uncertainty  in 
M.q,  is  ±10%  or  less.  Again,  as  with  altitude  for  a  given  flight  case,  the 


(97) 

(97a) 


temperature  as  a  function  of  altitude  would  be  known  to  much  better  than  ±10% 
and  the  resultant  uncertain/  in  MU, would  correspondingly  be  much  less. 


For  Mach  numbers  between  .  5  <  Mffi  <6,  it  was  seen  in  Section  II  that 
the  cosine  exponent  n  in  the  basic  relation  for  the  pressure  distribution  was 
Mach  number  dependent  (see  Figure  16).  In  addition,  it  was  shown  that  the 
value  of  Y\  could  be  calculated  once  ot  and  ^3  were  obtained.  Therefore, 
using  Figure  1  6  as  a  calibration,  it  would  be  possible  to  obtain  an  estimate  of 

in  this  transonic-supersonic  range.  It  is  felt,  however,  that  a  calibration 
better  defined  through  further  experiments  would  be  required  to  make  this 
Mach  number  determination  scheme  valuable. 


SECTION  y 


-’Sit 


ABLATION  AND  RESULTANT  UNCERTAINTIES 
ON  AIR  DATA  ATTITUDE  ANGLE  OUTPUTS 


1.  ABLATION  PERTURBATIONS 

Thus  far,  the  air  data  equations  have  been  derived  for  a  hemispherical 
probe.  However,  ablation  of  the  nose  can  result  in  a  change  in  nose  shape 
and  therefore  perturb  the  equations.  Given  the  re-entry  trajectory,  and  the 
nose  material  and  diameter,  the  amount  of  nose  recession  can  be  calculated. 
The  perturbation  on  the  equations,  due  to  ablation,  is  then  treated  as  an 
uncertainty  in  the  pressure  equation  exponent,  ±  AW ,  where  P=P5C.oS  0  . 
The  resultant  uncertainty  in  tX  or  y3  for  the  two  orifice  per  plane  calibration, 
Equations  (32)  and  (41),  is  then  determined  from  error  Equations  (36)  or  (42). 
While  by  no  means  rigorous,  this  analysis  was  performed  in  order  to  obtain 
an  order  of  magnitude  indication  of  the  possible  effects  caused  by  ablation. 

2.  ABLATION 

A  re-entry  trajectory  chosen  for  ablation  analysis  is  the  DODGO  trajectory 
presented  in  Reference  (20).  The  trajectory  was  extrapolated  down  to  60,  000 
ft.  tc  extend  the  ablation  calculations  over  a  greater  range  of  interest  (see 
Figure  38).  In  addition  to  the  trajectory,  other  assumed  inputs  were: 

1)  nose  diameter;  D  =  6.0  inches 
(felt  to  be  a  lower  limit  on  size) 

2)  nose  material;  "Graphitite  G" 

(a  material  for  which  experimental  ablation  data  taken  in  the  CAL 
Wave  Superheater  Supersonic  Tunnel  are  available) 

3)  zero  angle  of  attack  and  sideslip  over  the  entire  trajectory. 

The  results  of  the  ablation  calculations  are  presented  in  Figures  39  to  41. 

In  Figure  39,  the  nose  temperature  time  history  is  given.  In  Figure  40,  the 
stagnation  point  heat  transfer  rate  is  presented.  Both  of  these  curves  are 
fairing 8  through  point  by  point  integrations  (symbols)  that  were  hand  calculated 
for  the  trajectory.  The  nose  recession  A  X  (see  sketch)  as  a  function  of  altitude 
is  presented  in  Figure  41.  The  maximum  ablation  is  seen  to  be  approximately 
.  825  inches,  or  =.275. 


^/r 


3.  ABLATION  EFFECTS 


For  the  conditions  presented  in  Section  V,  2,  ,  it  was  seen  that  the  maximum 
value  of  stagnation  point  nose  recession  was  &X/  r  =  .  275.  The  following 
assumptions  and  methods  were  used  to  interpret  this  recession  distance  as  a 
perturbation  to  the  air  data  equations.  First,  as  previously  stated,  only  the 
case  at  zero  attitude  angle  was  considered.  Second,  it  was  assumed  that  as 
the  stagnation  point  recedes,  no  ablation  occurs  at  the  probe  shoulder  and, 
further,  that  the  nose  shape  can  be  approximated  by  an  ellipse,  see  sketch.  v 


For  an  ellipse,  the  coordinates  of  the  orifice  can  be  found  from 


where 


a  *  r-Al 
b  »  r 


Solving  Equation  (99)  for  X  and  taking  the  first  derivative  results  in, 

_  V  [  I  2  A  $  f  a $. 


Next,  it  is  necessary  to  assume  that  the  pressure  ratio  PA  at  X  ,  y 
can  be  approximated  by  a  cosine  function  even  though  the  surface  is  elliptical 
rather  than  circular. 


Thus, 


/ 


^/p  =  cos  (j)  -  cos"  0 


(101) 


where  </ V  is  the  angle  between  the  axis  and  a  line  normal  to  the  body  surface. 
Therefore, 

=  tarn  (f)  (102) 


Alternatively,  a  model  in  which  the  nose  shape  remains  spherical  but 
increases  in  radius  of  curvature  due  to  ablation  was  also  considered, 
however,  the  resultant  perturbation  to  the  air  data  equations  is  nearly 
the  same  as  the  above  analysis,  therefore,  it  is  not  included  herein. 


and  further 


(103) 


Substituting  Equations^  (102)  and  (103)  into  Equation  (100)  for  port  location, 
0  =  45  3  results  in  0  =  27°  45'. 


From  Equation  ( 1 C 1 ) 

Y\'  =  ■  ■  w  jf j -  (104) 

log  (cos  <p) 

and  since  the  h  value  used  in  the  assumed  distribution  is  ft  =  2.  24  (for  the 
hypersonic  case),  we  get  an  error  in  Y)  of: 

An  =  -  i. 44 

From  Section  II,  the  resulting  error  in  for  the  two  orifice  configuration  is, 


A  OL  - 


(coX  +‘taviotJ  — -  2C  +•  p °&e (cos 0  +  s • n 5^, ^ vs otjj i 


Values  of  Aot.versus  forAtt  =  1.  44  are  given  in  Figure  42.  Here,  values 
attx^Oare  considered.  To  arrive  at  these  conditions  physically,  it  would  be 
necessary  to  fly  at  CL=  0  while  the  nose  ablates  symmetrically,  and  then  go  to 
some  angle  of  attack  while  the  nose  is  still  symmetrically  ablated.  The  curve 
for  An  =  1.  44  shows  that  large  errors  exist  in  (X.  due  to  ablation.  However, 
this  perhaps  is  the  largest  possible  ablation  case  in  consideration,  since  the 
stagnation  point  was  kept  fixed  and  because  the  nose  diameter  considered  was 
the  lower  limit  in  size.  The  ablation  equations  (not  presented  here)  show  an 
ablation  dependence  on  nose  radius  of  the  form 


(105) 


Thus,  for  example,  for  a  1 2 -inch  diameter  probe 

» ,o<n 

r 

and  for  this  case,  Am  =  .  60.  The  angle  of  attack  error  A  oC  ,  for  D  =  12  inches 
is  also  shown  in  Figure  42.  While  a  12 -inch  diameter  probe  reduces  A  Ci 
significantly,  it  can  be  seen  that  the  perturbation  is  still  large.  Thus,  it  could 
be  concluded  that  nose  recession  should  be  limited,  either  by  cooling,  or  by 
large  nose  diameter  to  Ai/y  less  than  perhaps  5%  if  satisfactory  attitude  angle 
results  are  expected,  using  the  two-orifice  per  plane  model  and  equations. 

Alternatively,  the  sensitivity  of  the  attitude  equations  to  the  exponent  SI  can 
be  eliminated  by  using  three  pressure  ports  instead  of  two  in  each  attitude 
plane.  For  the  case  of  three  orifices  and  a  logarithmic  pressure  relation, 
Equations  (45)  and  (52),  it  is  possible  to  eliminate  the  pressure  distribution 
exponent  n  and  for  this  case  the  change  in  nose  shape  due  to  ablation  will  not 
cause  an  error  in  attitude  determination.  Pressure  and  attitude  equations  for 
this  orifice  configuration  are  described  in  Section  H.  It  should  be  cautioned 
however,  that  for  this  scheme  to  be  correct,  the  nose  pressure  distributions 
must  obey  a  cos  0  relation  even  though  the  value  of  V>  is  unknown. 
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SECTION  VI 


SUMMARY  OF  RESULTS  AND  RECOMMENDATIONS 


1.  SUMMARY  OF  RESULTS 


Air  data  outputs  obtainable  from  pressure  measurements  on  a  hemisphere 
probe  have  been  investigated  analytically  for  free  stream  Mach  numbers  which 
include  and  put  emphasis  on  the  hypersonic  flow  regime.  Specifically,  vehicle 
attitude  can  be  obtained  from  the  pressure  inputs  alone,  whereas  true  air  speed 
requires  the  additional  input  of  free  stream  density.  A  simple  flow  chart  of  air 
data  inputs  and  obtainable  air  data  outputs  at  hypersonic  and  supersonic  Mach 
numbers  are  shown  in  Tables  Ha  and  Hb,  respectively.  A  more  detailed  cnarfc, 
Table  HI,  lists  air  data  inputs,  outputs,  equations,  assumptions,  etc.,  and 
provides  in  effect  a  summary  of  the  study  performed.  At  hypersonic  Mach  ^ 
numbers,  vehicle  attitude  or  angle  of  attack  06  and  sideslip  can  be  obtained 
using  only  pressure  distributions  measured  on  a  hemispherical  nose  through 


the  relations 


(32) 

(41) 


Attitude  angle  uncertainties  in  these  expressions  resulting  from  one  percent 
individual  pressure  measurement  uncertainties,  and  possible  pressure  distribution 
variations  are  on  the  order  of  ±1/2  degree  for  attitude  angle  variations  from 
-20°  to  +50°,  The  equations  can  be  used  for  angles  as  high  as  85°  with  no 
resultant  discontinuities,  but  with  uncertainties  becoming  as  large  as  ±6°. 


Over  the  lower  Mach  number  regime  (transonic,  supersonic)  the  recommended 
Mach  number  independent  attitude  angle  expressions  are 


\og  (  Cos  ~  $ )/cos  oQ 
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(54) 


Attitude  angle  uncertainties  for  these  expressions  are  also  on  the  order  of 
±1/2°;  however,  the  useful  attitude  range  for  these  expressions  is  approximately 
±20°,  the  expressions  becoming  unusable  at  higher  angles  due  either  to  dis¬ 
continuities  or  multiple  value  solutions. 


Velocity  (true  air  speed)  can  be  obtained  for  Mach  numbers  of  approximately 
two  or  greater  from  the  hemisphere  pressures  if  free  stream  density  is  an 
additional  input  using  the  expression  ( 

1.02!  (%,) 
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(92) 
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Likewise,  indicated  air  speed  can  be  obtained  from  the  hemisphere  pressure 
inputs  using  y 

* ...  'rv 


u.  = 


(95) 


[cos  [W(fZ  (£)-M]  cos[(fr I) (^(Aj-F,))]  j 

where  /00  is  sea  level  atmospheric  density. 

For  Mach  numbers  i  (o}  K1  equals  2.  24  whereas  tt  is  Mach  number  dependent 
for  2<  M  <  (o. 


An  altitude,  or  more  properly,  a  pressure  altitude  cannot  be  obtained 
from  the  pressures  measured  on  a  hemisphere  in  hypersonic  flow.  Given  free 
stream  density,  a  density  altitude  can  be  obtained  to  within  at  least  ±8%. 

As  with  altitude,  the  pressures  measured  on  the  hemisphere  are  insufficient 
inputs  to  obtain  hypersonic  free  stream  Mach  number.  When  density  is  given, 
however,  and  velocity  is  calculable,  Mach  number  can  be  obtained  to  ±10%  or  better. 

For  the  simple  attitude  expressions,  Eq.  (32),  (41),  large  errors  in  attitude 
due  to  ablation  could  be  expected.  However,  use  of  Eq.  (47),  (54)  would  allow  attitude 
to  be  determined  even  with  ablation  if  the  nose  pressure  distribution  still  obeyed  a 
COS  n  0  relation. 

2.  RECOMMENDATIONS 


The  present  analysis  of  air  data  outputs  obtainable  •  t  hypersonic  Mach 
numbers  using  measurements  of  nose  pressure  has  led  ro  a  simplified  set  of 
air  data  equations.  These  expressions,  believed  to  hold  with  acceptable 
accuracy  in  the  high  Mach  number  flow  regime,  resulted  from  methods  used 
to  predict  the  hemispherical  pressure  distribution  over  a  broad  range  of  flight 
conditions.  Since  these  equations  allow  easy  utilization  of  the  air  data  inputs, 
they  appear  to  be  very  promising  for  hypersonic  flight  usage.  It  would  seem 
worthwhile,  therefore,  to  study  the  equations  experimentally.  That  is  not  to 
say  that  considerable  experimental  data  have  not  been  used  in  the  derived 
expressions;  however,  the  equations  have  not  been  experimentally  verified 
over  the  range  of  attitude  angles  (up  to  50®)  and  flight  velocities  (to  20,000 
feet  per  second)  considered  in  this  study. 


Specifically,  experimental  verification  of  the  two  following  expressions 
upon  which  the  air  data  equations  are  constructed  would  be  of  primary  interest. 


P  =  F,  cosn  &  n  -  2.  24 

ps  =  K  K  =  1.92 


(10) 

(106) 


These  tv/o  expressions  should  be  studied  experimentally  over  the  flight 
envelope  of  interest  with  variation  in  Uoo  ,  (altitude),  oL  and  .  In 

addition,  insofar  as  possible,  other  test  variables,  e.  g,  Mach  number,  Reynolds 
number,  model  size,  should  be  duplicated. 


Flight  duplication  of  velocity  presents  the  most  difficult  requirement  for 
ground  test  facilities.  With  regard  to  this  problem,  hypersonic  shock  tunnel 
facilities  appear  to  be  the  most  promising,  the  Cornell  Aeronautical  Laboratory 
Shock  Tunnels  for  example  having  the  capability  of  duplicating  velocity- dsn*- 'ty 
flight  conditions  for  a  full-size  probe  up  to  true  air  speeds  of  about  15,000  feet 
per  second  (see  Reference  21). 
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In  conclusion,  if  the  air  data  equations  for  OL  ,  /3  and  U®  can  be  shown 
to  hold  satisfactorily  over  the  flight  regime  predicted,  then  the  use  of  the 
hemisphere  probe  in  actual  vehicle  flights  is  recommended. 
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TABLE  31  b 


TRANSONIC  -  SUPERSONIC  OUTPUT  FLOW  CHART 
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